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Abstract 

Total Ordinary Muon Capture (OMC) rates are calculated on the basis of the Quasiparticle 
Random Phase Approximation for several spherical nuclei from ^''Zr to ^'^^Pb. It is shown 
that total OMC rates calculated with the free value of the axial-vector coupling constant gA 
agree well with the experimental data for medium-size nuclei and exceed considerably the 
experimental rates for heavy nuclei. The sensitivity of theoretical OMC rates to the nuclear 
residual interactions is discussed. 



1 Introduction 



The main aim of investigating muon capture on atomic nuclei is the determination of 
the couphng constants of the weak hadronic current in the nuclear environment. In 
principle one is interested in the weak axial- vector gA_ and pseudoscalar gp coupling 
constants. The observables of nuclear muon capture are calculated as functions of gA 
and gp and then compared to the corresponding experimental data. The Radiative 
Muon Capture (RMC) is traditionally considered as the most promising source of 
information on gp The sensitivity to gp could be increased if one considers the 
ratio of the total RMC rate to the total Ordinary Muon Capture (OMC) rate ^. 
Therefore the problem arises how to make consistent calculations of the total RMC 
and OMC rates within the same nuclear model. 

The investigation of total OMC rates is also of interest in itself because it allows 
the determination of g^ independent from beta-decay. In contrast to beta-decay all 
the final nuclear states with noticeable transition strength can be populated in muon 
capture. Therefore, small variations in the low-energy parts of theoretical strength 
functions give no dramatic changes in the calculated OMC rates in contrast to calcu- 
lations of log(/t) values. An additional advantage of muon capture as compared to 
beta decay is that the total OMC rates are measured for many stable and long-living 
nuclei. Therefore, one can in OMC not only study the variation of g^ with nuclear 
mass, but more delicate effects, such as the isotopic dependence of the effective weak 
interaction constants between leptons and nucleons. 

The present situation in the field of OMC rate measurements can be summarized 
as follows. Total OMC rates have been measured for many nuclei with high precision 
[^. The theoretical interpretation of the experimental information is, however, still 
controversial. The main reason lies in the necessity of a correct description of the 
nuclear response in both, OMC and RMC. 

The theoretical investigation of nuclear muon capture has a rather long history. 
Up to now, three different approaches have been developed for calculating total (ex- 
clusive) OMC rates on complex nuclei. The first approach is based on the closure 
approximation and related sum rule methods. In both cases the energy of the out- 
going neutrino is replaced by some average value which is a parameter of the theory. 
The OMC rate is then obtained using the closure relation for the final nuclear states 
The second approach utilizes the local density approximation. Here the OMC 
rates are calculated for infinite and uniform nuclear matter as a function of the proton 
and neutron densities. The OMC rate for a finite-size nucleus is then obtained by 
integrating this function over the realistic density distribution or by determining its 
value for a certain value of the nuclear matter density 0, ^ . The common drawback of 
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both approaches is that nuclear muon capture is considered without any connection to 
other processes which may occur in a nucleus. Also, the collective nature of the nuclear 
response to external fields is lost. The third approach, which is used in this paper, does 
not suffer from these defects. Here the (exclusive) OMC rates are calculated for all 
definite final states of the product nucleus and the total rate is obtained by summing 
over all considered final states. 

Atot = 5I%- (1) 
/ 

The calculations of Refs. 1^, 11, have been carried out within this approach. 

In the present work we study the total OMC rates on heavy nuclei within a mi- 
croscopic description of the nuclear structure. The one-body effective Hamiltonian of 
nuclear OMC obtained within the Morita and Fujii formalism jl^ is used. The wave 
function of the bound muon and its binding energy are calculated approximately [14|, 
taking into account the finite size of the nucleus. The nuclear matrix elements of the 
effective OMC Hamiltonian and the excitation energies of the states of a product nu- 
cleus are calculated within the Quasiparticle Random Phase Approximation (QRPA), 
an extension of the usual RPA to nonclosed shell nuclei. Velocity-dependent terms are 
included in the calculations and have been evaluated with Woods-Saxon single-particle 
wave functions having the correct asymptotic behaviour. 



2 The Effective Hamiltonian of Nuclear Muon Capture 

The total rate of OMC is calculated by summing the rates A ji for all partial transitions 
i ^ f. In a spherically symmetric nucleus Ajj is given by [p!5| ] 



Afi= [GcosQcYiE, 



Mi + m' 2Ji + 1 



^ " (2) 
Y.{Ml{u) + Ml{u + 1) + Ml{-u) + Ml{-u - 1)) 

u 

The M„(k)'s are the amplitudes for the transition in which a neutrino is created 
in a state with energy Ey and angular quantum number k {k = I for j = 1 — 1/2 
and K = —I — 1 for j = I + 1/2). u is the angular momentum transferred to the 
nucleus. These amplitudes are combinations of the weak form factors with nuclear 
matrix elements 
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Mu{u) = J^^(^^Gv[Ouu]-^I^^Ga[Iuu] 
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Mu{-u-l) = J^1-^(^./^[^Gv[0uu] + ^Ga[Iuu] 



-^-[lu + lup] 



u + l 



2u + 

Here the effective form factors are defined in the usual way 

Gv = 9v{q^){l+Ej2M)+gs{q^) 

Ga = gA{q^)-{9v{q')+9M{q^)){Eu/2M) (4) 
Gp = {gp{q^)-gA{q^)-9T{q^)-9v{q^)-9M{q^)){E,/2M) 

and the nuclear single-particle matrix elements are given by 

A 

[OH = {JfWsJT-Y. M^k) ttju{E,rk) Y^ifk) II Ji)/^'^Jf + 1 

fe=l 



k=l 

= i{Jf II sj -7-^^^{rk)tl jyj{E^rk) \y'w{fk) ®Pk]u \\ Ji)/^2Jf + 1 

fc=i 

[O'U'up] = i{Jf II Y 2~ X] ^/^(''fc) ju{E^rk) Yuifk) {Pk,^k) \\ Ji)/ \J'^Jf + '^ 

k=i 



(5) 



where = v'7r/2xJ„+i/2(2;) is a spherical Bessel function and 

[cr(8)y^(f)]u,„i„ = ^ (lmu;mt„|um„)c7^yt„,„i^(f) (6) 



is the tensor product of two spherical tensor operators. i?!>^(r) is the radial wave function 
of the bound muon. For the nucleon isospin operators t+ the convention t'^\p) = \n) 
is used. 
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Eqs.(|2|) and (|5|) show that the capture rate depends rather strongly on the energy 
of the outgoing neutrino 

F ( I I ^ rn^-\eis\+ Mj- Mf - E* 
= {m, - \e,s\ +M.-Mf-E)(l ^(m, + M.) ^ 

For a large Z nucleus the muon binding energy ei^ has to be calculated taking into 
account the finite size of the nuclear charge distribution. The excitation energy of the 
final nuclear state E* enters into the capture rate mainly through the neutrino energy. 
In order to obtain the transition energies and nuclear matrix elements one has to 
use some nuclear model. In the present work the QRPA has been used since it gives 
a precise prescription of how the interaction between particle and hole excitation can 
be included in the calculation of the nuclear matrix elements and the excitation 
spectrum of the product nucleus. 



3 The Nuclear Model 

A detailed description of the nuclear model used in the present work can be found in 
and will not be repeated here. The nuclear Hamiltonian consists of a mean field 
part, a monopole pairing interaction between like particles and a residual interaction. 



Hm — ^ ( -f^mean(i3) + -f^pair(i3) ] + -f^resid (8) 
i3=±l/2^ ^ 



For protons and neutrons separate Woods-Saxon potentials including spin-orbit inter- 
actions have been used to approximate the mean field. The mean field and pairing 
Hamiltonians 

j,m jm, j'm' , , 

(9) 

are in a first step approximately diagonalized by the special Bogoliubov transformation 
(see e.g. 

Solving the BCS equations leads to the Independent Quasiparticle Hamiltonian 

^^0(^3) ^ ^ eit3«imt3"i™*3 (11) 

jm 

with 

G 



^jt3 = V (^J*3 - + Cl and Ct^ = ^Y. ^jt^v^t^ ■ (12) 



5 



An averaging over the quasiparticle vacuum state ajmtslO) = is implied. As residual 
interaction we use an effective isospin-invariant separable interaction of the form 



I E (4'' + 4''(n-r2))QL^,(l)QijM(2) . 



^ L,J,M 

where Qlm and Qljm are single-particle multipole and spin-multipole operators and 
the isospin structure is displayed explicitly. Using 

(ri • T2) = 4(fi • t2) = 4i°t° + 2{t{4 + tft^). (14) 

one can combine the isospin operators with the multipole and spin-multipole operators 
leading to a new set of single-particle operators 

QiAiif^) = E {fT^'t'S^ fL{rk)YLM{k)tl\jmh)a\,^,^,^ajrnu (15) 

and 

QljmW= E {fmt'^\i^fLj{rk)[YL{k)(g)a{k)]jMt'j^\jmt3)a],^,^,ajrnt3 (16) 

3 

where is an isospin operator out of the set {i^, = t^, = tf± it\, (k = 1, 2)}. 
In the literature one finds two most frequently used variants of the radial form 
factors /L(r) and /Lj(r). These are 

fLir) = f{r)Lj = r'^. (17) 

and 

hir) = f{r)LJ = f{r) = ^U{r), (18) 

where U{r) is the central part of the shell-model potential used in Hmean- The mixed 
products tit2 and tit2 in Eq.(14) lead to particle-hole excitations changing the charge 
of the nucleus and are therefore involved in a description of charge-exchange interac- 
tion processes such as /3-decay, /x-capture or direct {p,n) and (n,p) reactions. The 
corresponding parts in the Hamiltonian (^) are constructed with the singl-particle 
operators 

^JM= E {jnmn\OjMt^\jpmp)a]^„^^aj^^mp (19) 

and their Hermitian conjugates, where Ojm can be either i'^ fj{r)YjM{r) or 
fLj{r)[YL{r) ®(t]jm respectively. 

The residual interaction ( [T^ ) contains only scalar products of the form 
{[Yj-i{fi) ® [Yj-iih) ® (T2]j) and {[Yj+i{fi) ® [^7+1(^2) ® ct2]j). The 
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tensor interaction which would mix the [Yj^i{fi)0ai]jM and [yj+i(r2) (8><T2] jm terms 
has been omitted because its inclusion would slightly influence the properties of the 
charge-exchange resonances in the range where the coupling constants assume reason- 
able values. 

To achieve an approximate diagonalization of the residual interaction Hamiltonian 
in the QRPA one introduces phonon creation and destruction operators. They are 
defined as linear combinations of tensor products of the quasiparticle operators ajmts 
and ctj'm't' and their Hermitian conjug ates (Eq.|lD. 

JpiJn 

The phonon amphtudes ^j^j^ and (p'j^jp are orthonormahzed according to 

^(mO- E{^;..C-CCJ = ^m'- (21) 

Having expressed the residual interaction in terms of the phonon operators ( pO| ) one 
obtains the QRPA equations through a variational principle pBI 



sf^mMHM^W/l) - {\HM\)-u;,mi,i) - 1)} = (22) 

using the normalization of the phonon amplitudes as a subsidiary condition. |) is the 
phonon vacuum: ^ij^l) = 0, approximating the ground state of a double even nucleus. 



The QRPA equations resulting from (22) are a system of homogeneous linear equations 
determining the phonon amplitudes tpj^j^^ and </'j„jp and the excitation energies iOi. 



-ujig\ + Rlq,w\, =0, 
where the following abbreviations have been used 



(23) 



5g = V'i,j„+Cin' «^^ = Cjn-Cjn' <l^ijp^3n), (24) 
eg = ^iv + ^in , = ± Vj^Uj^ (25) 

Kq' = ^iKq' - i^J^ (^i^ihiu^ K^'t' + i^i^hf^u^ K'^^i) (26) 
for naturaly parity states and 

for unnatural parity states, and h^^ stand for the reduced matrix elements of the 



multipole (15) and spin-multipole ( [Iq) single-particle operators with = t . The 
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QRPA amplitudes for the transitions from the even-even ground state to the excited 
states with total spin JM, and energies uji ± (A„ — Ap) are given by 

1 



Jp-3n 



if the charge of the nucleus decreases from Z to Z — 1 (as, for example, in {n,p) 
reactions) or by 



... 1 



if the charge increases from Z to Z + 1 during the transition (as in {p,n) reactions). 

In the following we will make a few remarks on the nuclear model parameters. 
There are three kinds of parameters which characterize the nuclear structure of our 
model Hamiltonian. The parameters of the mean field (Woods-Saxon potential) are de- 
termined such that they reproduce best the single-particle excitations of neighbouring 
odd-mass nuclei. Different potentials have been used for protons and neutrons. The 
parameters of the monopole pairing interaction have been obtained from the odd-even 
effect of nuclear binding energies. Since the theoretical binding energies depend on the 
single-particle energies there is a certain correlation between the pairing constants and 
the parameters of the Woods-Saxon potential. Usually one set of monopole pairing 
constants and Woods-Saxon potential parameters has been used in the calculations for 



a whole group of neighbouring nuclei. These parameter sets have been taken from [17|. 
The residual interaction between particle-hole excitations causes the collective small 
amplitude vibrations. The effective constants of this interaction have been determined 
by comparing the calculated excitation energies and the transition strengths of the 
collective isovector states with experiment. In the present work only the isovector in- 
teraction is of interest since we are dealing with a charge changing process. A detailed 
discussion of the used form of separable residual interaction can be found in the work 
of [^]. In order to demonstrate the sensitivity of the calculated total muon capture 
rate Atot to the residual interaction, both variants of radial form factors have been 
used with several values for the isovector coupling constant. 
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4 Results of the Calculations 



In this section we present detailed results of the calculated total OMC rates obtained 
for spherical nuclei of different mass regions. The data are presented in Table || to 
Table |l^. The capture rates shown in Fig. |l| and Fig. ^ are presented in the form of 
running sums 

A{E)= ^k, (30) 

k■.E^,<E 

where the energies are measured with respect to the ground state of the initial nucleus. 
In order to get a feeling for the sensitivity of the theoretical muon capture rates on 
the chosen shape and strength of the residual interaction, all calculations have been 
performed for both types of interaction ( p^ ) and (p^). Throughout we have employed 
gp/dA = 6.0 in the calculations. 



4.1 ""Zr and "^Mo 

The distribution of the strength of the transition operator at^ over the excitation en- 
ergies (Gamow- Teller strength function) has been studied in detail for the reactions 
9°Zr(p, n)9°Nb [|l9|, ^ and ^^Zt{^U, ^Ref^Nh |2|. The strength function has a promi- 
nent peak at an energy of 15.6 MeV. For this peak the reduced probability of the GT 
transition B[GT) = 10. The total observed transition strength below 20 MeV excita- 
tion energy is around 20. The total GT strength should be larger than 3(A^ — Z) = 30, 
the value given by the Ikeda sum rule p^ . 

Table |^ shows how the calculated properties of the GT strength function depend on 
and nf-. The model of non-interacting quasiparticles (residual interaction switched 
off: = = 0) cannot give the correct position of the peak of the strength 
function. For L = and J = 1, the residual interaction ( |T7| ) reduces to the simple 
(cj, (?) interaction considered in [23|. The calculations with = —23 /A recommended 



in correctly reproduce the position of the maximum of the strength function. 
The results of calculations with slightly different constants = —25/A and = 
—28/^4), presented in Table ||, show that the position of the maximum of the strength 
function is sensitive to k^^. The strength in the peak region does not depend on k^^ 
and considerably exceeds the experimentally observed one. So, the best value of the 
effective constant for the residual interaction (p7|) is = — 23/yl. 

Table || shows that the calculated energy of the GT resonance is less sensitive to 
for the residual interaction with form factor (|l8|). The correct position of the 
resonance is reproduced with = —0.33/A. Simultaneously the strength below 
and in the resonance region is considerably smaller than the corresponding strength 



calculated with interaction (17) (24.6 to be compared to29.2) and is closer to the 
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experimental value. The rest of the transition strength is located at high-excited l"*" 
states [p^ ]. Recently, new experimental data of the GT transition strength at very 



high excitation energies became available |2C]. The total B{GT) strength calculated 
with the residual interaction (fsl) is equal to 32.0. This agrees with the experimental 



value 34.2 ± 1.6 obtained in [20| by a multipole decomposition of the experimental 
^^7ji(j), n)^''Nb cross sections. From this consideration one concludes that the residual 
interaction (^) provides a better description of the at~ strength function than the 
interaction (0). 

Recent measurements of the af^ transition strength in ^''Zr(n,p)^'^Y [25| can also 
be compared with our calculations, since in this reaction the charge of the nucleus 
decreases as in muon capture. The crt^ strength, summed over all experimentally 
observed states is B^{GT) = 1.0 it 0.3 [25|. Our calculation with the interaction 
( |l7|) gives the following distribution of the transition strength. A considerable part 
of B^{GT) is concentrated on the first 1+ state of ^°Y. The other l"*" states have 
excitation energies between 10 and 15 MeV. For each of these states, B~^{GT) < 0.2. 

The af^ strength function calculated with the interaction (|l^) differs from the 
strength function obtained with interaction (|T7|). The B^{GT) calculated with in- 
teraction (|T8[) is almost three times the value of B^iGT) obtained with interaction 
(p!7|). This is due to highly excited states which are absent in the calculations with the 
/(r) = interaction. As before, the strongest transition goes to the first 1^ state of 
^'^Y, but the strengths of transitions going to 1^ states with energies between 5 and 
15 MeV are comparable to B^{GT). In this case the transition strength is distributed 
more uniformly over the excitation energies and the shape of the theoretical strength 
function is closer to the experimental one. 

It should be noted that the energy of the first l"*" state calculated for each of the two 
residual interactions does not depend on , and the corresponding Bf{GT) decreases 
slightly whith growing \ki^\. This indicates that already in ^°Zr the neutron excess 
prevents the creation of low-lying collective 1^ states in ^"Y. The residual interaction 
with /(r) = dU/dr ( |l8|) couples single-particle states with wave functions having the 
same orbital quantum numbers and different number of nodes in the radial parts. Due 
to the residual interaction (p!8|), these particle-hole excitations interact with particle- 
hole states which are all members of one spin-orbital multiplet and create in this way 
high-excited collective states |24]. The transitions to those states increase B^{GT). 

It was shown in that the theoretical OMC rates are rather insensitive to the 
constants of the multipole residual interaction k/. Nevertheless, the constant of the 
isovector monopole residual interaction, can be determined from analyzing isobar 
analog states (IAS). Results for the O"*" charge-exchange excitation in ^''Zr are presented 
in Table |2|. The independent quasiparticle model describes the 0^ charge-exchange 
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states as a set of non-interacting two-quasiparticles states. The states carrying the 
main transition strength are (059/2)^(059/2)11 and (lc?i/2)p(l^^i/2)n for the t~ or (p, n) 
branch and {0fY/2)p{0f 1/2)11 and (0d5/2)p(l(i5/2)n for the t+ or {n,p) transitions. Both 
interactions produce a collective state in the (p, n) excitation branch whose energy 
coincides with the experimental IAS energy of 12.0 MeV @ (k? = -0A3/A for inter- 
action (|l^ ) and = —35.0/A for interaction (|T7|)). In both cases this state consumes 
almost all of the t~ transition strength. 

In the {n,p) branch the interaction ( [TtI ) with the above mentioned coupling con- 
stant is not capable to produce a collective state and the strength goes into the 
(0/7/2)^(0/7/2)71 and (0(i5/2)p(ld5/2)n two-quasiparticle states. The interaction (|T^) 
with the above mentioned coupling constant is strong enough to produce a collective 
state at approximately the energy of the two-quasiparticle states. The difference be- 
tween the total t~ and transition strengths is constant and does not depend on the 
residual interaction. Therefore one can use the strength functions to determine the 
isovector residual interaction coupling constants and k^. For L > we have used 
the relations 

^01 

Ki'^ = "^"iry interaction (17) (31) 

and 

n^-J = for the interaction (18) (32) 

In our theoretical total OMC rates final states with J'^ = 0^, 1^, 2^ and 3^ have 
been taken into account. Contributions from final states with J > 3 turned out to 
be less than 1%. In Table ^ and Fig. || we present the calculated total OMC rates for 
^^Zr. The rates shown in Fig. |^ were calculated with both types of interaction, using 

= -0.43/^ and = -0.33/^ for type and = -35/A and = -23/A 
for type (p!7|). It can be seen that the main difference between the rates for the two 
different residual interactions originates from muon capture populating highly excited 
states. This difference is biggest for the 1"*" final states. Both calculated total rates 
agree with each other to within less than 5%. However, larger differences appear in the 
partial rates feeding states with a specific J'" . In general, one observes that for ^4 ~ 90 
nuclei the total capture rates calculated with the interaction (|l^ is larger than that 
calculated with the interaction (|T8[). 

The experimental total capture rate of A = 86.6.10^s~^ [Q] is in rather good agree- 
ment with both theoretical values. However, one should keep in mind that the ex- 
periment has been done with the natural isotope composition of ^^Zr, whereas the 
calulations refer to a specific isotope. 

The theoretical capture rates Atot show only a very slight sensitivity to a variation 
of gp/ga- The variation in A^t for ^^Zv amounts to 85.5.10^s~^ to 76.8.10^s~^ if 
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Qpl QA is increased from 4.0 to 12.0. The contribution of the velocity dependent matrix 
elements [l^ynp] and [Onnp] to Atot is rather small. The change in Atot is less than 2% 



when they are switched off. The estimated contribution of (12 - 15)% made in [11| 
could thus not be confirmed by the direct calculation. 

The interesting observation was made that the l"*" state around 5 MeV which gives 
a prominent contribution in the ^''Zr(n,p)^''Y reaction, contributes only very little to 
the total muon capture rate. This is pointing to the fact that one cannot directly 
use matrix elements from (n,p) or (p,n) reactions to analyze muon capture data. A 
large B(GT) does not always imply a large capture rate Ag. Here the explanation is 
that the radial matrix elements with spherical Bessel function jo{ki^r), dominating in 
allowed 0"*" — > l"*" partial transitions, are suppressed by the centrifugal barrier of the 
two-quasi particle state {0g7/2)n{0gg/2)p- This different appearance of the l"*" state of 
the product nuleus in various reactions finds its explanation only within a microscopic 
nuclear model. 

Table |^ and Table ^ contain the results of our calculation for ^^Mo. There is a 
satisfactory agreement between the experimental data and the rates obtained with our 
microscopic model. As in the case of ^^Zv the two different residual interactions lead to 
a different population of the specific states which after summing over all final states 
is smoothed out. A comparison with the measured capture rate on natural Mo seems 
reasonable, since ^^Mo is the lightest even Mo-isotope, and the capture rate decreases 
with increasing neutron excess. 

4.2 OMC on Even Tin Isotopes 

In this section we consider the even tin isotopes ^^^~^^^Sn, a long chain of stable spheri- 
cal nuclei. Tin isotopes have a completely filled OgQ/2 proton subshell and gradually 
filled OgT/2, lc^5/2) ^d-^/2 2si/2 neutron subshells. Table ^ shows the total OMC 



rates calculated with two sets of single-particle potential parameters |17]. It can be 
seen that the dependence of the total capture rates on the Woods-Saxon potential 
parameters is stronger than its dependence on the residual interaction coupling con- 
stants. The calculated rates show also a strong dependence on the neutron excess of 
the target nucleus. 

The experimental value of the capture rate measured for natural Sn is 106.7 • 10^ s"'^ 
[Q]. Considering that ^^^'^^'^Sn contributes more than 50% to the natural isotope mix- 
ture, shows that the Woods-Saxon model parameters fitted to ^ = 121, Z = 51 give 
a better description of the total rates. We can say that the agreement between the 
theoretical capture rates and the experiment is quite good. 
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4.3 Heavy Nuclei with Large Neutron Excess. ^^°Ce and ^°*Pb 

The most important observation to be made in the case of heavier nuclei is (see Fig. ^ 
the large theoretical overestimation of the rates. The difference between the calcu- 
lated total OMC rates becomes larger with increasing mass number. Table |^ shows 
the results for ^^'^Ce and Table ^ those for ^'^^Pb. The rates calculated with the resid- 
ual interaction (ITtI ) are smaller than those obtained with the interaction (|T8|). This 
difference comes about mainly due to capture populating the high-excited l"*" states 
which is absent in calculations using the interaction (^^. The experimental energies 
of the collective IAS and GT {at~) states can be reproduced with the interaction ( |l7|) 
using = -28.0/^ to obtain 18.94 MeV for the IAS state and using nf = -23.0/^ 
to obtain the GT state at an energy of 19.71 MeV. The corresponding experimental 
energies are 18.8 MeV for the IAS state and 19.2 MeV for the GT state respectively. 
More than 80% of the total GT strength sits at the peak of the distribution. Using the 
interaction (18) one obtains the collective 1^ state at 16.85 MeV using n^^ = — 0.43/yl. 
Approximately 50% of the total GT strength goes to this state and more than 30% of 
the strength is shifted to the higher l"*" states. These 1^ states with high excitation 
energies are responsible for the fact that the collective 1^ state remains in a region 
below 18 MeV even if \tif\ is doubled. 

The experimental value of htot for ^osp]-, jg j^gg . j^qS g-i g_ Thus both calculations 
overestimate the total rate considerably. It is therefore interesting to compare our 
results with those from previous work achieving a good agreement with experiment 

|ll|. There is, however, a difficulty in performing this comparison because the 
authors of [0, 10, 11] presented their results as relative contributions to A^t, related 



to a specific angular momentum transfer L. To compare these data to our results, 
we have to assume that the transitions to the 0^ and 1^ final states proceed via the 
L = transition. The transitions to the 0^, 1^, 2^ states are accompanied by a L = 1 
orbital angular momentum, etc. With this assumption we implicitely assume that 
|[1J - 1J]| >> |[1J + 1J]|. A test shows that deleting of [IJ + IJ] reduces the 
corresponding Ajj by less than 10%. The fractional contributions obtained by this 
procedure are given in Table ^. A comparison with the results of 10, shows 
that the main difference to our results comes from the contributions of the 0^ — > 1^ 
transition. This prominent role of the 0^ 1^ transition was also found in a recent 
calculation |]l^]. The total OMC rate for ^'^^Pb obtained in the present work with the 
interaction ( p!7| ) compares well with the value Atot(^'^^Pb) = 161 • lO^s^^ obtained in 
[ p^ ] using a 5-function residual interaction. 
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5 Discussion and Conclusion 



In this work a theoretical evaluation of the total OMC rates for medium - weight 
and heavy spherical nuclei using QRPA was presented. For the first time an attempt 
was made to include the velocity-dependent terms, evaluated with single particle wave 
functions having the correct asymptotic behaviour. It was shown that the contribution 
of these terms to Aj^^Jr™ is rather small. To avoid confusion, some remarks about the 
meaning of "velocity-dependent" terms have to be made. Usually all terms having its 
origin in the small components of the nucleon 4-spinors are called "velocity-dependent" 
terms. However, in the derivation of the effective muon capture Hamiltonian [^, part 
of these terms experience a transformation due to momentum conservation 

p + ^ = n + v, 0, =^ n = p — P. (33) 

As a result, only the gradient acting on the proton wave functions is left. We have 
explicitly shown that the matrix elements with the proton gradient give a minor con- 
tribution to the capture rate. So one can conclude that the main effect of the velocity- 
dependent terms is already accounted for in the effective Hamiltonian due to momen- 
tum conservation. 

Our calculations show that the total OMC rates are not very sensitive to the con- 
stants of the nuclear residual interactions. On the other hand they may strongly depend 
on the shape of the residual interaction used in the calculations. The main influence 
on Atot calculated with different residual interactions comes from the difference in the 
description of GT transitions, O'*' ^ 1^, especially at higher excitation energies. 

A comparison of our theoretical total OMC rates with experiments shows (Table |lO| 
and Fig. |2|) the following situation. For the medium weight nuclei (^'^Zr, ^^^~^^^Sn) a 
reasonable agreement between theory and experiment can be achieved using the free 
values of gA and gp. No renormalization of g^ is needed in this mass region. 

The AJJqj"'' exceed, however, considerably the experimental values for the heavier 
nuclei ^^'^Ce and ^'^^Pb. Therefore, in order to reproduce the experimental values, a 
renormalization of gA seems to be necessary for heavy nuclei. This renormalization 
is model dependent; it depends in our case on the coupling constants and the shape 
of the residual interaction, as can be seen from Table |^ and From Fig. |2| one is 
tempted to deduce some systematic deviation of Aj:^™'' from the experimental data. 
This makes it impossible to draw any definite conclusion on the necessity of a quenching 
of gA for heavier nuclei. It seems that the nuclear model used in this investigation 
reaches its limits of useful application so that further theoretical studies are necessary. 



The widespread belief |11] that a theoretical description of total OMC rates faces no 



particular problem seems to be untenable. 
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Table 1: Properties of charge-exchange l"*" excitations in ^^Zr calculated with two variants of residual interactions. 





at 


as in 


[p, n) reaction 


crt+ as in 


(n,p) reaction 


Energy of 

maximum 


B-{GT) 


Energy of 

maximum 


s+(Gr) 


Icjlal 


ill max. 


V)(;knv max. 


lot al 


ill max. 


l")('l(nv max. 


0.00 


11.84 


32.06 


16.96 


13.79 


5.34 


3.32 


2.21 


0.00 










/(r) = dU/dr 








-0.23 


14.93 


31.89 


20.88 


6.23 


5.62 


3.15 


1.31 


0.00 


-0.33 


15.76 


32.03 


19.88 


4.67 


5.68 


3.30 


1.14 


0.00 


-0.43 


16.36 


32.21 


18.41 


3.74 


5.73 


3.47 


1.03 


0.00 










fUr) = r 


L 








-23.0 


15.73 


30.63 


23.83 


5.32 


5.63 


1.89 


1.12 


0.00 


-25.0 


16.09 


30.56 


23.98 


4.86 


5.65 


1.82 


1.07 


0.00 


-28.0 


16.63 


30.05 


23.86 


4.28 


5.67 


1.72 


0.99 


0.00 



Table 2: Properties of charge-exchange 0"'" excitations in ^^Zr calculated with two variants of residual interactions. 





t- 




Energy of 
maximum 


B- 


Energy of 
maximum 


B+ 


total 


in max. 


below max. 


total 


in max. 


below max. 


0.00 


4.62 


10.14 


8.62 


0.00 


16.17 


0.29 


0.07 


0.13 










/(r) = dU/dr 








-0.33 


10.99 


10.58 


8.30 


0.37 


18.08 


0.72 


0.55 


0.16 


-0.43 


12.00 


10.70 


7.94 


0.14 


18.34 


0.84 


0.67 


0.16 


-0.53 


12.81 


10.81 


7.44 


0.07 


18.56 


0.96 


0.79 


0.15 










fUr) = r' 










-23.00 


10.35 


10.04 


7.87 


1.78 


16.91 


0.18 


0.05 


0.07 












16.20 




0.04 




-25.00 


10.71 


10.03 


8.29 


1.32 


16.92 


0.18 


0.05 


0.07 












16.20 




0.04 




-28.00 


11.28 


10.02 


8.68 


0.88 


16.92 


0.17 


0.05 


0.07 












16.20 




0.04 




-31.00 


11.87 


10.02 


8.88 


0.62 


16.92 


0.16 


0.05 


0.06 












16.20 




0.04 




-34.00 


12.46 


10.01 


8.96 


0.45 


16.93 


0.15 


0.05 


0.06 












16.20 




0.04 




-37.00 


13.06 


10.00 


8.97 


0.34 


16.93 


0.15 


0.04 


0.05 












16.21 




0.03 





Table 3: The rates of OMC (in 10^ s ^) on ^^Zr summed over the final states with 
specific spin and parity J'^ . The second line gives the contribution of each final state 
in%. 





final states, J'" 


total 
rate 


0+ 


0" 


1+ 


1- 


2+ 


2" 


3+ 


3^ 


0.00 


4.9 


3.6 


24.3 


43.4 


15.2 


19.3 


11.6 


3.3 


125.6 




3.9 


2.8 


19.3 


34.6 


12.1 


15.4 


9.2 


2.6 


100.0% 










f{r) 


= dU/dr 








-0.23 


5.3 


2.2 


28.3 


27.3 


8.5 


12.3 


5.5 


1.4 


90.8 




5.9 


2.5 


31.2 


30.0 


9.3 


13.5 


6.0 


1.6 


100.0% 


-0.33 


5.3 


2.2 


29.1 


25.1 


7.4 


11.2 


4.4 


1.2 


85.8 




6.2 


2.5 


33.9 


29.2 


8.6 


13.1 


5.2 


1.4 


100.0% 


-0.43 


5.3 


2.2 


29.9 


23.5 


6.6 


10.5 


3.7 


1.0 


82.8 




6.4 


2.6 


36.1 


28.4 


8.0 


12.7 


4.5 


1.2 


100.0% 




















-23.0 


4.7 


1.9 


23.4 


27.2 


10.2 


12.2 


7.1 


1.9 


88.7 




5.3 


2.2 


26.4 


30.6 


11.5 


13.8 


8.0 


2.2 


100.0% 


-25.0 


4.7 


1.9 


23.2 


27.2 


10.0 


11.8 


6.9 


1.9 


86.9 




5.5 


2.2 


26.7 


26.5 


11.5 


13.6 


7.9 


2.2 


100.0% 


-28.0 


4.7 


1.8 


23.0 


25.7 


9.6 


11.3 


6.5 


1.8 


84.4 




5.6 


2.1 


27.2 


30.4 


11.4 


13.4 


7.8 


2.1 


100.0% 
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Table 4: Properties of charge-exchange l"*" excitations in ^^Mo calculated with two different nuclear residual 



interactions. 





at 


as in 


[p, n) reaction 




at' 


^ as in 


(n,p) reaction 




Energy of 
maximum 


B-{GT) 


Energy of 
maximum 


B+{GT) 


total 


in max. 


below max. 


total 


in max. 


below 


max. 


0.00 


11.85 


28.79 


16.94 


10.70 




4:.U 


5.96 


4.75 




0.00 










/(r) = dU/dr 












-0.23 


14.98 


27.42 


11.61 


12.17 




5.04 


4.59 


2.92 




0.00 


-0.33 


15.61 


27.28 


16.94 


4.83 




5.19 


4.45 


2.57 




0.00 


-0.43 


16.15 


27.23 


16.50 


3.42 




5.29 


4.40 


2.33 




0.00 










fL{r) = r 


L 












-23.0 


15.54 


26.36 


18.36 


6.86 




5.07 


3.53 


2.56 




0.00 


-25.0 


15.81 


26.23 


19.40 


5.61 




5.11 


3.40 


2.45 




0.00 


-28.0 


16.25 


26.06 


20.20 


4.47 




5.16 


2.23 


2.29 




0.00 



Table 5: The rates of OMC (in 10^ s ^) on ^^Mo summed over the final states with 
specific spin and parity J'^ . The second line gives the contribution of each final state 
in%. 





final state J'" 


total 
rate 


0+ 


0" 


1+ 


1- 


2+ 


2" 


3+ 


3^ 


0.00 


5.3 


4.1 


27.4 


52.2 


18.1 


24.2 


13.5 


3.8 


148.7 




3.6 


2.8 


18.4 


35.1 


12.2 


16.3 


9.1 


2.6 


100.0% 










f{r) 


= dU/dr 








-0.23 


5.6 


2.5 


30.5 


33.3 


10.1 


15.6 


6.5 


1.7 


105.8 




5.3 


2.4 


28.8 


31.5 


9.5 


14.7 


6.2 


1.6 


100.0% 


-0.33 


5.6 


2.4 


31.0 


30.6 


8.8 


14.1 


5.3 


1.4 


99.2 




5.6 


2.4 


31.3 


30.8 


8.8 


14.2 


5.4 


1.4 


100.0% 


-0.43 


5.6 


2.4 


31.5 


28.7 


7.8 


13.2 


4.5 


1.2 


95.0 




5.9 


2.6 


33.2 


30.3 


8.3 


13.9 


4.8 


1.8 


100.0% 




















-23.0 


5.1 


2.3 


26.2 


33.6 


12.3 


15.9 


8.5 


2.3 


106.2 




4.8 


2.1 


24.6 


31.6 


11.6 


15.0 


8.0 


2.1 


100.0% 


-25.0 


5.1 


2.2 


26.0 


32.9 


12.0 


15.5 


8.2 


2.2 


104.1 




4.9 


2.1 


24.9 


31.6 


11.5 


14.9 


7.9 


2.2 


100.0% 


-28.0 


5.1 


2.1 


25.7 


31.8 


11.6 


14.9 


7.9 


2.2 


101.2 




5.1 


2.1 


25.4 


31.5 


11.5 


14.7 


7.8 


2.1 


100.0% 
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Table 6: Total OMC capture rates on Sn isotopes (in 10^ s ^). 



Target 




SW parameters 




SW parameters 


nucleus 


for (ID 


115,49 


121,51 


for (ID 


115,49 


121,51 




-0.23 


139.2 


130.1 


-23.0 


141.7 


123.0 




-0.33 


130.5 


123.1 


-25.0 


138.9 


120.8 




-0.43 


124.9 


119.0 


-28.0 


135.2 


117.4 


ii^Sn 


-0.23 


130.0 


122.1 


-23.0 


131.9 


115.0 




-0.33 


122.1 


116.1 


-25.0 


129.4 


112.7 




-0.43 


117.1 


112.4 


-28.0 


125.9 


109.5 


i20Sn 


-0.23 


121.2 


111.8 


-23.0 


122.6 


107.3 




-0.33 


114.2 


109.5 


-25.0 


120.3 


104.1 




-0.43 


109.8 


106.3 


-28.0 


117.2 


102.1 


i22Sn 


-0.23 


113.2 


107.8 


-23.0 


114.0 


99.7 




-0.33 


106.9 


103.2 


-25.0 


118.0 


97.7 




-0.43 


103.0 


100.6 


-28.0 


108.0 


95.0 


i24Sn 


-0.23 


105.5 


101.3 


-23.0 


105.5 


91.7 




-0.33 


99.9 


97.1 


-25.0 


103.4 


89.9 




-0.43 


96.5 


95.1 


-28.0 


100.7 


88.3 
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Table 7: The rates of OMC (in 10^ s -'^) on ^'^'^Cc summed over final states with specific 
spin and parity . The second line gives the contribution of each final state in %. 





final state 


total 
rate 


0+ 


0- 


1+ 


1" 


2+ 


2- 


3+ 


3" 










f{r)-- 


= dUjd 


r 








-0.23 


11.7 


3.8 


55.4 


39.9 


20.2 


21.0 


10.3 


3.5 


165.7 




7.0 


2.3 


33.4 


24.1 


12.2 


12.7 


6.2 


2.1 


100.0% 


-0.33 


11.7 


3.5 


57.4 


38.0 


17.9 


21.3 


8.7 


2.9 


161.2 




7.2 


2.2 


35.6 


23.5 


11.1 


13.2 


5.4 


1.8 


100.0% 


-0.43 


11.7 


3.4 


59.2 


36.7 


16.2 


21.7 


7.7 


2.5 


159.0 




7.3 


2.1 


37.3 


23.1 


10.2 


13.7 


4.8 


1.6 


100.0% 










fdr) = 










-23.0 


9.8 


3.4 


41.4 


35.8 


25.2 


15.1 


12.9 


5.1 


148.9 




6.6 


2.3 


27.8 


24.1 


16.9 


10.1 


8.7 


3.4 


100.0% 


-25.0 


9.8 


3.3 


41.2 


35.1 


24.7 


14.8 


12.5 


5.0 


146.3 




6.7 


2.2 


28.1 


24.0 


16.9 


10.1 


8.6 


3.4 


100.0% 


-28.0 


9.8 


3.1 


40.8 


34.1 


23.9 


14.4 


12.0 


4.8 


142.8 




6.9 


2.2 


28.5 


23.9 


16.7 


10.1 


8.4 


3.4 


100.0% 
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Table 8: The rates of OMC (in 10^ s^^) on ^O^Pb summed over final states with specific 
spin and parity J'^ . The second line gives the contribution of each final state in %. 





final state 


total 
rate 


0+ 


0- 


1+ 


1" 


2+ 


2- 


3+ 


3" 










f{r)-- 


= dUjd 


r 








-0.23 


15.7 


3.1 


70.5 


33.8 


29.7 


29.5 


10.4 


5.8 


198.6 




7.9 


1.6 


35.5 


17.0 


15.0 


14.9 


5.2 


2.9 


100.0% 


-0.33 


15.7 


2.7 


72.5 


33.1 


26.9 


31.8 


9.5 


4.9 


197.2 




8.0 


1.4 


36.8 


16.8 


13.7 


16.1 


4.8 


2.5 


100.0% 


-0.43 


15.7 


2.4 


71.9 


32.6 


25.0 


33.2 


9.0 


4.3 


194.2 




8.1 


1.2 


37.0 


16.8 


12.9 


17.1 


4.6 


2.2 


100.0% 










fdr) = 










-23.0 


13.5 


2.8 


46.9 


23.8 


33.7 


19.9 


11.0 


8.0 


159.5 




8.5 


1.7 


29.4 


14.9 


21.1 


12.4 


6.9 


5.0 


100.0% 


-25.0 


13.5 


2.6 


46.7 


23.3 


33.0 


19.8 


10.7 


7.8 


157.4 




8.6 


1.7 


29.7 


14.8 


21.0 


12.6 


6.8 


5.0 


100.0% 


-28.0 


13.5 


2.4 


46.5 


22.6 


32.0 


19.6 


10.3 


7.5 


154.5 




8.7 


1.6 


30.1 


14.7 


20.7 


12.7 


6.7 


4.9 


100.0% 
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Table 9: Fractional contributions of different multipoles to Atot for ^osp]-, (gj^her con- 
tributions of the orbital momentum transfer or contributions of the transitions to the 
states with specific J'"). 



Reference 


AL 











1 


2 


3 




i 




26 




14 




48 


12 


1 


10 


1 


29 




13 




52 


7 


1 


y 


1 


23 




34 




34 


8 


(a) 






45 




~ 35 






18 


~ 2 


(b) 




38 




~ 29 






27 


~ 5 
















0+ 


1+ 


0" 


1" 


2" 


2+ 


3+ 


3" 


1 


n 


1 


5 


36 


6 


14 


12 


28 


4 


2 


(a) 




8 


37 


2 


17 


16 


14 


5 


3 


(b) 


8 


30 


2 


15 


12 


21 


7 


5 



(a) - present paper, calculations with the residual interaction (|T8|), ^f"^ • A = —0.43; 

(b) - present paper, calculations with the residual interaction (|T7|), • A = —25.0. 
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Table 10: Summary of results : total OMC rates (in 10^ s ^) calculated for gp/gA = 
6.0 with two different radial form factors of the nuclear residual interaction. 



Target 
nucleus 




A for dU/dr 


K^-^A for 


expr. 

i 


-0.23 


-0.33 


-0.43 


-23.0 


-25.0 


-28.0 


90Zr 


90.8 


85.8 


82.8 


88.7 


86.8 


84.4 


86.6 


92Mo 


105.8 


99.2 


95.0 


106.2 


104.1 


101.2 


92.2 




130.1 


123.2 


119.0 


123.2 


120.8 


117.4 




iiSSn 


122.1 


116.1 


112.4 


115.0 


112.7 


109.5 




i20Sn 


114.8 


109.5 


106.3 


107.3 


105.1 


102.1 


106.7 


i22Sn 


107.8 


103.2 


100.6 


99.7 


97.7 


95.0 




i24Sn 


101.3 


97.1 


95.1 


91.7 


89.9 


88.3 




i40Ce 


165.7 


161.2 


159.0 


148.9 


146.3 


142.8 


114.4 


208pb 


198.6 


197.2 


194.2 


159.5 


157.4 


154.5 


134.5 
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w 
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Figure 1: ^"Zr : Total and partial integrated capture rates up to the excitation energy 
E of the final nucleus. The partial rates are shown for the final states with J'^ = 1^. 
Solid lines: results of the calculation with the residual interaction ([T^), dashed lines: 
the same with the interaction (^) , the solid vertical bars show the distribution of the 
calulated partial rates over the excited states with J'^ = 1^. In the upper part of the 
figure the distribution over all excited states is shown. 
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Figure 2: Comparison between experimental and calculated total capture rates. 
Crosses: experimental data taken from Q and Q (the vertical extension indicates 
the experimental error), full circles: calculation with the interaction (18), open circles: 
calculation with the interaction (|l^) . The full curve connects the rates calculated with 
the Goulard - Primakoff formula O with its parameters determined from M. 
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Figure 3: 208p]-, . Total and partial integrated capture rates up to the excitation energy 
E of the final nucleus. The partial rates are shown for the final states with J'^ = 1^. 
Solid lines: results of the calculation with the residual interaction ([T^), dashed lines: 
the same with the interaction (^) , the solid vertical bars show the distribution of the 
calulated partial rates over the excited states with J'" = V^. In the upper part of the 
figure the distribution over all excited states is shown. 
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